DUDLEY KNOT HDD. ART 
NAVAL PC j’iT"' i , • ■ i.’F SCHOOL 
MONTEP.Pv • . ..:v ? • 



:ooc* 



NAVAL POSTGRADUATE SCHOOL 

> 

Monterey, California 




THESIS 

SENSITIVITY AND ROBUSTNESS ANALYSIS 
FOR SEA-SKIMMING MISSILE 

by 

Olavo Amorim de Andrade 
September 1986 

Thesis Advisor: D. J. Collins 

Approved for public release; distribution is unlimited. 



T 23 00 33 




UNCLASSIFIED 

SECURITY CLASSiEfCAfiQN QF~fmS PAG? 



REPORT DOCUMENTATION PAGE 



la REPORT SECURITY CLASSIFICATION 

UNCLASSIFIED 



lb. RESTRICTIVE MARKINGS 



2a SECURITY CLASSIFICATION AUTHORITY 



2b. (^CLASSIFICATION /OOWNGRAOING SCHEDULE 



3 DISTRIBUTION/ AVAILABILITY OF REPORT 

Approved -for public release; 
distribution is unlimited. 



4 PERFORMING ORGANIZATION REPORT NUMBER(S) 



S. MONITORING ORGANIZATION REPORT NUMBER(S) 



6a. NAME OF PERFORMING ORGANIZATION 

Naval Postgraduate School 



6b. OFFICE SYMBOL 
(If applicable) 

67 



7a. NAME OF MONITORING ORGANIZATION 

Naval Postgraduate School 



6c. AODRESS (City, State, and ZIP Code) 

Monterey, Cali-fornia 93943-5100 



7b. AODRESS (City, State, and ZIP Code) 

Monterey, Cali-fornia 93943-5100 



Ba NAME OF FUNDING / SPONSORING 
ORGANIZATION 



Bb. OFFICE SYMBOL 
(If applicable) 



9. PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 



8c. AO OR ESS (City, State, and ZIP Code) 



10 SOURCE OF FUNDING NUMBERS 



PROGRAM 


PROJECT 


TASK 


WORK UNIT 


ELEMENT NO. 


NO 


NO 


ACCESSION NO 



1 1 TITLE (include Security Classification) 

SENSITIVITY AND ROBUSTNESS ANALYSIS FOR SEA-SKIMMING MISSILE 



12 PERSONAL AUTHOR(S) 



Andrade, Olavo A. 



13a TYPE OF REPORT 


13b TIME COVEREO 


14 0ATE OF REPORT (Year, Month, Day) 


15 PAGE COUNT 


Engineer ' s Thesis 


FROM TO 


1986 September 


120 



16 SUPPLEMENTARY NOTATION 



17 



COSATI CODES 



18 SUBJECT TERMS ( Continue on reverse if necessary and identify by block number) 



FIEL0 


GROUP 


SUB-GROUP 


Sensitivity, Eigenvalues, 


Robustness , 








Aerodynamic Coe-f f icients, 


Sea-Skimming Missile, 








Optimization, Minimum Singular Value 



9 ABSTRACT ( Continue on reverse if necessary and identify by block number) 

This thesis presents the application o-f EIGENVALUE SENSITIVITY 
ANALYSIS and SINGULAR VALUE ANALYSIS to the control o-f a SEA-SKIMMING 
supersonic missile, in the vertical plane. 

The study is divided in -four basic parts: 

a) The development o-f the model; 

b) EIGENVALUE SENSITIVITY ANALYSIS with respect to the variation o-f 
the aerodynamic parameters o-f the autopilot/air-frame o-f the 
mi ssi 1 e; 

c) Analysis o-f the time response with respect to the 
variation o-f the aerodynamic parameters; and 

d) Robustness analysis and improvement o-f the system, using the 
SINGULAR VALUE ANALYSIS. 



20 OtSTRlBUTlON/ AVAILABILITY OF ABSTRACT 

XlUNCLASSIFlED/UNLIMITEO □ SAME AS RPT. □ OTIC USERS 



21. ABSTRACT SECURITY CLASSIFICATION 



Daniel J. Collins 



22 b TELEPHONE (Include Area Code) 

(408) 646-2826 



22c OFFICE SYMBOL 
67Co 



ID FORM 1473. 84MAR 



83 APR edition may be used until exhausted 
All other editions are obsolete 
1 



SECURITY CLASSIFICATION OF THIS PAGE 



UNCLASSIFIED 



UNCLASSIFIED 

SECURITY CLASSIFICATION of This PAGE Dmt m Enfrmd) 



19. ABSTRACT Continued 



All the analysis 
programs using the 
POSTGRADUATE SCHOOL. 



is based in results of simulation 
software available at the NAVAL 



S/N 0102- LF- 014- 6601 



2 UNCLASSIFIED 



SECURITY CLASSIFICATION OF THIS PAGEf»7i*n D»tm Entmrmd) 



Approved for public release; distribution unlimited. 



Sensitivity and Robustness Analysis 
for Sea-Ski mming Missile 



by 



□lavo Amorim de Andrade 
Li eutenant-Commander , Braz i 1 i an Navy 
B.S., Escola Naval , R. J ., Brasi 1 , 1970 
B.S. , Universidade de S=(o Paulo, S.P., Brasil, 

Submitted in partial fulfillment of the 
requirements for the degree of 

AERONAUTICAL ENGINEER 



from the 



1976 



NAVAL POSTGRADUATE SCHOOL 
September 1986 



ABSTRACT 



This thesis presents the application of EIGENVALUE 
SENSITIVITY ANALYSIS and SINGULAR VALUE ANALYSIS to the 
control of a SEA-SK IMM ING supersonic missile, in the vertical 
plane. 

The study is divided in four basic parts: 

a) The development of the model. 
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variation of the aerodynamic parameters of the 
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I. INTRODUCTION 



The majority of tactical missiles used against surface 
targets are SEA-SKIMMING and fly as close as possible to the 
surface of the water in order to make difficult the 
detection and reaction of the enemy. 

This work addresses the problem of analysing the design 
of the altitude control system of a supersonic SEA-SKIMMING 
missile in both aspects of SENSITIVITY and ROBUSTNESS, with 
respect to the variations of the aerodynamics parameters. 

Those variations can develop from uncertainties in the 
model as well as from changes in the flight conditions of 
the vehicle. 

The SENSITIVITY analysis makes it possible to verify, 
for a given model which parameters are more important for 
the desired flight path. 

The ROBUSTNESS analysis has the objective of verifying 
if the system could be affected by perturbati ons or noise 
and will be the base on calculations of modified feedback 
gains which are capable of producing the desired flight of 
the missile despite perturbat ions. 

The contents of the different sections are as following: 

Section II presents the model and the state equations of 
the system as well as the calculations to complete the chain 
of control of the missile, considering the restrictions of 
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the SEA-SKIMMING scenario, where the altitude control system 
cannot tolerate any overshoot to avoid that the missile hits 
the water. 

Section III develops the SENSITIVITY ANALYSIS taking 
into consideration the variation of the eigenvalues and time 

response of the system of state equations for variations of 

the aerodynamic parameters in a range of +257. of the 

original value. 

In section IV is described the ROBUSTNESS ANALYSIS based 
upon the minimum singular value of the return difference 
matrix as function of the frequency; section V introduces a 
design technique using optimization routine in order to 
obtain a "complete robust" solution and section VI presents 
the final conclusions. 
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ALTITUDE(FEET) 



1 1 . MATHEMATICAL MODEL OF THE MISSILE 



The work assumes a SEA-SKIMMING missile with the -flight 
in the vertical plane according to the path shown in Figure 
II. 1. 



O 

O 




Figure II. 1 - Missile Trajectory in the Vertical plane 

The missile is assumed to fly at constant velocity, at 
150 ft, when it is commanded to a final altitude of 30 ft. 
Only the pitch channel will be analysed. 
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A. THE MISSILE 

The missile will be considered as presented by Arrow 
CRef.ll. This model has been selected in order to avoid 
classification problems. The geometry as well as its 
characteristics are presented in Appendix A. 

The states that were considered in our model are defined 
in Table I as well as the relevant aerodynamics parameters. 

In the Figure II. 2 we have the schematic representat i on 
of the uncoupled pitch control chain of the missile, 
neglecting the yaw and roll movements; 



H 

s 




Figure II. 2 - Block Diagram of the Pitch Control Chain of 
the Missile. 



where 

T) , — "bias" accel erat i on ; 

zcb 

T» - commanded accel erati on ; 
zc 

1 )^ - achieved acceleration; 
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v(t> - velocity; 

H(t) - altitude; 

H — commanded altitude, 

s 

The Transfer Function G(s) represents the missile itself 
and the other blocks are the controller. The input to this 
block will be the "commanded 11 acceleration and the output 
will be the achieved accel eration. 

The commanded acceleration comes from the feedback of 
the states of the controller (compensated state, velocity 
and altitude) added to a "bias" accel erat i on . 

In the diagram of Figure II. 3, we have a complete 

representati on of the system with the states used in 
modelling our missile. 

The "STATE EQUATIONS" which represent our system are the 
foil owi ng : 



x = -150.0 x 4 - 2,646.66 x - 705.75 x^_ 

11 5 fr 

x = 1 . 353 x + x , - 1 . 353 u + 1 . 334 H 

2 13 ZC S 

x = -6.572 x - 5 x, + 6.572 u -6.4B1 H 

3 13 zc s 

x = -44.332 x — 59.109 x 

■ 56 

X = x u - 0.1482 x — 0.0395 x, (II. 1) 

5 “ 5 o 

x fc = -188.4 x 7 + 188.4 x # 
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x 7 = -0.4608 x t - 2.231 x £ - 0.3406 x^ + 2.231 x h + 

-15.0949 x _ 20.13 x, - 0.143 x, + 0.4608 t) 

s b 7 'zc 

- 0.4544 H 

s 

x = -50.0 x o - 495.0 x + 495.0 H 

8 8 9 S 

x 4 „= -17.644 x - 4.705 x 

10 S b 



This system of equations can be represented in matri 
■form, where we have: 



x = A x + B u 



A 



-150 0 0 

1.35 0 1.0 

-6.5 0 -5.0 

0 0 0 

0 0 0 

0 0 0 

-0.46 -2.2 -3.4 
0 0 0 

0 0 0 

0 0 0 



0 

0 

0 

0 



1.0 



0 



2.2 



0 

0 

0 



2.646.66 -705.74 O O 



0 

O 

-44.33 
-O. 15 
O 

-15. 1 
0 
0 

-17.64 



O 

O 

-59. 1 1 
-0.04 
-188.4 
-20. 13 
0 
0 

-4.71 



0 0 
0 0 
O O 
O 0 
188.4 O 
-0. 14 O 
O -50. 
0 O 
0 O 



0 O 
0 O 
O O 
O 0 
O O 
O 0 
0 0 
-495 0 
0 1 
0 0 
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B 



O O 

- 1.353 1.334 

6.572 - 6.481 

O O 

O O 



O O 

O O 

0.461 - 0.454 

0 495.0 

0 0 

u = — F x 



F = 



0 0 0 0 0 0 0 
0 0 0 0 0 0 0 



0.049 0.98615 1.404 
0 0 0 



The vector input "u" has two components 

i) - the commanded acceleration; and 
zc ’ 

H - the desired altitude. 

5 

The input i> z( _ is considered as two parts; one 

corresponds to the -feedback o-f the states (from the 
controller) and another is the “bias" acceleration of lg. 
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TABLE I 

DESCRIPTION OF THE SYSTEM STATES AND 
AERODYNAMIC COEFFICIENTS 



1 o 



& 



pc 



filtered acceleration 

acceleration compensator 

acceleration compensator 

q - pitch angular rate 

a — angle of attack 

pitch tail incidence 

commanded pitch tail incidence 
... x - controller compensator 

H(t) - altitude 

v(t) - velocity 



AERODYNAMIC COEFFICIENTS 



C 

C 

C 

C 



fTloC 



slope of curve of pitching 



, change in C per degree pitch 
m&p m 



No: 



slope of curve of normal force 



change in per degree pitch 



moment coef f i c i ent 

control incidence 

coefficient C x 
N 

control incidence 
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U.<t> , 1 v ( t ) , 1 H(t) 




Figure II. 3 - The Missile Model 
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(oas/pej) b 



The model can be separated in two main parts; the 
CONTROLLER and the AUTOPILOT /AIRFRAME. 

B. THE CONTROLLER 

The controller was developed based in the procedure 
introduced by Dowdle [Ref. 23. 

The missile was considered as a point mass and, 
therefore the altitude as function of time is calculated by 
double integration of the achieved accel erat i on . 

The controller was developed without take in 
consideration the missile itself, i.e., the transfer 
function G(s) in Figure II. 2 was assumed as "unity". 

According to the [Ref. 23, this procedure is valid if the 
bandwidth of G(S) is larger than the crossover frequency of 
the controller. 

Due to the constraints of the scenario, the controller 
has to be designed in such way that no overshoot will be 
tolerate. Further requirements are that the angl e-of -attack 
and the accel erati on , have to remain below a certain level. 

The controller are responsible for the commanded 
acceleration by feedback of its states. 

The technique used to solve the feedback problem was the 
LINEAR QUADRATIC STATE FEEDBACK REGULATOR. This solution has 
a large stability margin with respect to phase variations, 
as well infinite gain margin. 
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The "cost function" to be minimized is the followings 



J = f < zlt) Q z(t) + ult, R u(t> 3- dt (II. 2) 
^ 0 

In Figure 1 1. 4 we have the schematic representat i on of 
the CONTROLLER. 



D 



z 



H 

s 



1 


V(t) 


i 


H ( t ) 

+ ' 


1 

1 V 


G (s) 
c 


s 




s 


■V_y 



Figure 1 1.4 - The Controller 



The LEAD 
in order to 
Its transfer 

G (s) 
c 



compensator indicated 

improve the transient 

function is s 

K*(s+w ) 
z 

(s+w ) 
p 



in the controller is used 
response characteristics. 



where K = 10.0 

w = 0.5 rad/sec 
z 

w = 50.0 rad/sec 
P 

In order to have the state equations of the CONTROLLER, 
we draw in FIGURE II. 5 its signal flow graph using the 
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adequate transformation of the compensator and showing the 



states that will be feedback. 



f 




Figure II. 5 Signal Flow Graph of the Controller 

The state vector is defined as: 

x ( t ) = C x(t) (H(t)-H ) v ( t ) 1 T (II. 2) 

c s 

In the Figure II. 5 the states of the controller are 
represented as indicated in Table I. 

From the diagram, we have the following STATE EQUATIONS: 

x = A x + B u (II. 3) 



z = D x 
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The state matrices are: 



A = 



w K ( w — w ) O 

P z P 

0 0 1 

0 0 0 



B = 



O 

0 

1 



D = C 1 K O 3 



Since we consider 
equation I I . 1 



to feedback all the states, let in 



Q = I (Identity) 



and 



R = 0.04 



Using the OPTSYS program to solve that OPTIMUM CONTROL 
problem, the following feedback gains were obtained; 

F l = 0.049 

F 2 = 0.98615 

F 3 = 1.4044 

correspondents to the feedback of the states indicated 
in Figure II. 5. 
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C. THE AUTOPILOT /AIRFRAME 



The controller was developed assuming the transfer 
■function o-f the Autopi 1 ot /Air-frame as unity. This assumption 
implies that the achieved acceleration is equal to the 
commanded acceleration; however, they are different due to 
the aerodynamic of the missile. In order to take that in 
consideration we have to introduce the autopilot/airf rame 
model . 

Normally, the sea-skimming is a SKID-TO-TURN missile 
(turns instantaneously), but due to problems of 
classification we have used for our Sea — Skimming model the 
Autopilot given by Arrows CRef.ll for a BANK-TO — TURN 
missile (has to bank in order to turn) that was designed to 
fly at 30,000.0 ft instead of sea level. 

Take in consideration that the main objective of this 
work is to establish a procedure for sensitivity and 
robustness analysis and that we are studying only the 
vertical movement of the missile uncoupled of any other 
movement, we consider that the approximation is useful. 

In Figure II. 6 we have the block diagram of the missile, 
representing the plant without the controller. 

The state equations derived from this diagram are on the 
Appendix A. 
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Control Law Actuator Dynamic Model 
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D. THE TIME RESPONSE 

Using the software CONTROLS, the system represented by 
equations (II. 1) was simulated and the time response for the 
altitude, velocity, acceleration and angle of attack were 
obtai ned . 

The altitude is in Figure II. 7 and we can see that the 
missile arrives at the desired altitude in about 10 seconds. 

No overshoot is shown and a small steady-state error is 
presented characterizing our Type 0 system. The steady-state 
error was less than 5/C. 

The angl e-of -attack , reproduced in Figure II. 8, varies 
between 3.7 deg and —1.2 deg satifying the condition 
required by the airframe, as indicated in CRefs. 1 and 23. 

The acceleration is in Figure II. 9 and the velocity is 
in Figure II. 10 representing a smooth behavior of the 
system. 

In the next chapter will be presented the Sensitivity 
Analysis using the developed model. 
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ORIGINAL SYSTEM 




(L-Oacmunv 



Figure II. 7 - A1 1 i tude <-f t ) vs. Time(sec) 
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4.0 6.0 8.0 10.0 12.0 14.0 16.0 

TIME(SEC) 



ORIGINAL SYSTEM 



o 




(oaa)XDvnv-io-aiDNV 



Figure II. 8 - Angle-of-Attack (deg) vs. Time(sec) 
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ORIGINAL SYSTEM 



o 




(OSS/l^AlIDCnSA 



Figure II. 10 - Velocity (-ft/sec) vs. Time(sec) 
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TIM1(SEC) 



III. SENSITIVITY ANALYSIS 



A. INTRODUCTION 

The objective of the SENSITIVITY ANALYSIS is to 
analyse how the behavior of the system is affected by 
variations of the aerodynamic parameters. Those variations 
are due to uncertainty of the design and/or problems ocurred 
during the flight. 

From the developed model we can see that the important 

parameters to our system are C . C.. , C and C., , as 
K 7 m<*’ No:' m & N£’ 

described on TABLE I. 

The variations could be interpreted as perturbat i ons in 
our STATE EQUATION SYSTEM (A-matrix > and the analysis, 
basically, take into consideration the change in the poles 
positions due to change in the parameters. 

As we have feedback, we replace the A matrix by the 
augmented A-matrix ( A_^ = A - B F ). 

The sensitivity of a matrix where the elements are 
ti me-invar i ant , can be verified in different ways. Among 
those procedures, some are enumerated below. 

First, using the method introduced by Golub and Van Loan 
CRef.31 we can analyse the matrix to verify how it is "ill- 

conditioned"; where an "ill-conditioned" matrix means a very 
sensitive system. 
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In order to accomplish the analysis, we have to 
calculate the "condition number". It is defined as: 

K( *V = 11 1 1 ii a -1 1 1 

o a a 

where 

I I | | is the norm of augmented A and 

I I A || is the norm of inverse augmented A 

Q 

The condition mumber depends upon the considered norm 
but it can be proven that the condition of the matrix will 
be the same for all norms. 

When we use the 2-norm one has 

K (A ) = * (A a } (III. 1) 

o (A ) 

— a 

where a (A ) and o (A ) are the maximum and minimum 
a - a 

singular value of the A-augmented matrix, respecti vel y . 

Calculating the condition number for the different 
values of the parameters it is possible to establish how the 
sensitivity of the system changes when there are changes in 
some elements of the state matrix. 

The problem of this kind of analysis is that we can not 
verify the sensitivity of the system with respect to a 
specified element of the state matrix but only the 
sensitivity of the overall system. 

In Frank CRef.43, we have another method where the 
analysis is done based upon the STATE SENSITIVITY EQUATIONS 
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that could be characterized in time domain, -frequency domain 



orin terms of a performance index. 

The basic idea is to define "sensitivity functions" and 
analyse the behavior of the system by those functions. 

The analysis starts with the system defined by: 



where : 

- x^ represents the state vector; 

- t is time ; 

- u is the input vector ; and 

- at represents a vactor of the different parameters. 

If at is changed by Ao: (some perturbati on > and using 
Taylor's expansion, we have: 



The subscript <* o indicates that the partial derivative 

is taken at the nominal parameter values. 

The partial derivative is called the trajectory vector ( 

x ) with respect to the j — parameter and has the same 

dimension as the state vector. 

The components of that vector are defined as the 
TRAJECTORY SENSITIVITY FUNCTION. 



x = f (>£, t , u , <*) 



(III. 2) 




J 



(III. 3) 



X. . (t,<* ) 

1 j ’ - o 



dx ( t , at) 

dot . 

J 



o 



(III. 4) 
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The equation III. 4 represents the partial derivative of 

the i state with respect to the j^- parameter. Assuming we 
have n states and m parameters, it will be possible to form 
a n x m matrix called TRAJECTORY SENSITIVITY MATRIX (III. 5). 
The matrix is also called the JACOBIAN MATRIX. 



X 



Jx_i 
dot , 



Jot 

m 



Jxn 

Jot 

1 



Jxn 

Jot 



m 



(III. 5) 



The change in the behavior of the system will be 
determined by: 

r 

Ax = X (t,oc ) Aot = y X. Aot . (III. 6) 

-0 - Z_, J J 

j = i 

The analysis is based in how the sensitivity function 
changes when we change the parameters. The result could be 
visualized by plotes of X as function of time, as states of 
the expanded system. 

Finally, we consider the EIGENVALUE SENSITIVITY 
presented by Porter & Crossley CRef.53. 

This latter procedure was choosen to be used in the 
analysis of this thesis, where the sensitivity of the sistem 
is related to the variation of the eigenvalues of the 
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augmented matrix (A = A - B F) with respect to the 

di 

variation of the parameters. 

In order to further develop our analysis, the concept of 
MODAL ANALYSIS will be introduced. 



B. THE MODAL ANALYSIS 

The concept is based in generating the input vector of a 
system by linear feedback of the state vector in such way 
that the prescribed eigenvalues are associated with the 
dynamic model of the resulting close-loop system. 

As an example, let us assume a scalar system: 

x = a x (t) + b u(t) ; (III. 7) 



if u ( t ) = 0 ====> x = x, exp (at) 

0 

where exp (at) defines the "mode" of the system. 

If a ^ 0 the system is stable and if a ^ 0 the system is 
unstable. 

Assuming g as the feedback gain, we have: 

u(t) = g x(t) (III. 8) 

Inserting the value of u from equation III. 8 in equation 
III. 7: 



x = (a + b g) x (t) (III. 9) 

and, solving for x(t): 

x(t) = x (0) expC (a + bg)tl 

Now the mode of the system is defined by the exp(a+bg). 
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When we have a state vector , the system is represented 
by 

x = A x + B u (III. 10) 

Considering 

u = - F x 

and assuming all the states as observable and 
control 1 ab 1 e 7 the substituting value of u in equation 
(III. 10) , we have 

x = (A - B F) x or x = A x (III. 11) 

a 

The precise nature of the free motion of the continuous 

time system following any disturbance can be described in 

terms of eigenvalues and eigenvectors of the augmented plant 

matrix (A ) . 

a 

Assuming A^ has n distinct eigenvalues ( X ± , X 2 , . . . , X r ) ; 

then it also has n corresponding linearly independent 
eigenvectors u lf u fif ... f u such that 

A u. = x. u. : i=l,2,...,n (III. 12) 

a l l l 9 77 7 

Using the calculated eigenvalues and eigenvectors we can 

write the modal matrix of A as: 

a 

U = Cu u . . . u ] 

12 n 

where each column of U is the eigenvectors of A^. 

Based in the modal matrix U, the equation (III. 12) could 
be written as 
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A U 
a 



A U 



(III. 13) 



with A as a diagonal matrix of the eigenvalues 



The matrix A is obtained -from equation (III. 13) 



U 1 A U 
a 



A 



(III- 14) 



In addition to the eigenproperties of A , the 

o 

correspond! ng properties of its transposed matrix (A') play 

a 

an important role in the modal analysis. 

The matrix A^ has the same eigenvalues of A^, but a 

different sets of eigenvectors.lt can be represented as 




J = 1,2 



1 ■ ■ ■ 1 n 




In CRef.5D it is shown that 



v . u . = u . v 



U . U . V 

J 1 1 J 




(III. 16) 



where &. . is the Kronecker delta: i.e. 

* j ’ 




if i = j 



and 




if i * j 



Thus in matrix form 



V'U = I 



(III. 17) 



where the collumns of V are the vectors v. 

i 



Using the modal matrix of A^ and introducing a new state 



vector y(t) , we have 



x (t) = U y (t) 



(III. 18) 
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The derivative of both sides, gives us: 

x (t) = U y ( t ) (III. 19) 

with x(t) = A x(t) (III. 20) 

Combining equations ( 1 1 1 . 18) , ( 1 1 1 . 19) and (III. 20), we 
have: 

y = if 1 A a U y (t ) 
and using the equation (III. 14), 



y (t ) = A y (t ) (III. 21) 

Because the matrix A is diagonal the equations are 
uncoupled and the solutions o-f the system represented by 
(III. 16) are given by 

y. (t) = y. (0) exp(X.t) (III. 22) 

7 1 7 l ^ l 

Aplying equation III. 18, 



x . (t) = [u, u_ ... u 3 
l 12 n 



(III. 23) 



which implies that 

x(t) = u^y 1 (0) exp(Xj't) + u 2 y 2 (0) exp(X 2 t) + 



Taking the value of y(t) in equation (III. 18), 

y (t) = U _1 x (t) (III. 24) 

Using equation (III. 17), 

y . (0) = v.' x (o) (III. 25) 

7 1 l 
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Combining equations (1 1 1.23) and (II I. 25): 



x (t) 




x(0) u . v ' exp(X.t) 



V . CA U \ A . 

11 ^ 1 



(III. 26) 



Equation 1 1 1 . 26 shows that the system represented by 
(III.Q) has its modes described by the eigenvectors of 



order to have the system stable, x(t) > 0 as t > ® , 

which happens if and only if 



In the Modal Control CRef.53 all this procedure is shown 
for the case when some of the eigenvalues are equal which 
does not apply to our system. 

As the role of the eigenvalues is so important, we can 
verify how the sensitivity of our control 1 system by 
analysing the change in the eigenvalues of the matrix A or 
the augmented A for the close loop system when any element 
is changed. 

In order to avoid calculating the new eigenvalues for 
each different A, the equation (III. 26) can be used to find 
approximated values that can be analysed analytically. 



and its transpose and eigenvalues of A^ matrix. Also 



Re X. < 0 

l 



for all the values of i 
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C. EIGENVALUE SENSITIVITY 

In our system we consider the change in the elements of 
the matrix due to variations in the aerodynamic 
coe-f -f icients. 

Basically, we concerned with the eigenvalues that are 
changed in the direction of the right hand plane, making the 
system unstable or affecting the time response. 

The result of the analysis could be plots or tables 
relating the eigenvalues with the variation of the 
parameters. 

For each change in the parameters, we have a new matrix. 

Approximations to the eigenvalues associated with the 
different plant matrices in the vicinity of the original one 
could be used without having to compute the eigenvalues of 
each new matrix. 

We describe the approach presented in the reference 
CRef.5D applying the "first order eigenvalue sensitivity" to 
the change of one of the parameters and doing the complete 
analysis with the calculated eigenvalues. 

1 . First-order Eigenvalue Sensitivit y 



From equation III. 11, we have 



x (t) = A x (t) 

Q 



Let 's assume 




where a. . is any element of A 
kl a 
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Differentiating equation (III. 12) with respect to 



a, , indicates that 
kl 



aA u. + A dm = axi + x aui 

aa. _ aa. . aa, . aa. . 

kl kl kl kl 



( 1 1 1 . 27 



multiplying (III. 27) by then gives 
v' aA u. + v A aui = v.' axi u. + x. v.' aui 

l—i l l - — l li 

ia ki aa ki aa ki aa ki 



As v.' A is equal to X. v.' , we have 

l l l ’ 



v.' aA u . = v.' axi u. 

l—i i l 

aa ki aa ki 



Which reduces to the set of equations 
axi = v. (k) u.(l) ; i=k=l=l ,2, . . . ,n (III. 29) 

— — l l ’ ’ 

* a kl 
Due to 



aA = a . , a . . 
al lk Jl 

aa ki 



In equation (III. 28), v^(k) and u^(l) represent the 

k— element of v^ and the 1— element of u^ , respectively. 

Those coefficients may be considered as the 
elements of a set of n eigenvalue sensitivity matrices. 

The new eigenvalues (x^) of the altered matrix, 

when the element a is changed will be calculated as 



A i v 

X . = X . + _i A a. _ 

l l t kl 

aa ki 
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When we have more than one parameter altered, the 
linearity could be used and the -final eigenvalue will 
obtained -from the summation of the sensitivity coefficients 
multiplied by the respective change in the coefficients. 



As an example of this procedure, let's consider our 

system, changing the value of C such that 

mjp 

AC, =0.02 
m&p 

When C , is changed as specified, the following 
m hp 

elements of the augmented-A matrix of the state equations 

presented in the Chapter II are changed: 

A (4,6): Aa u = 14.778 

a ’ h 6 



A (7,6): Aa_, = 5.035 
a 7 6 

The eigenvalues and eigenvectors of the augmented A 
matrix and of its tranpose were calculated using the program 
Mat lab. 



A 

a) Eigenvalue : 

X ± = -159.75 + j 18.85 

A X =Cv (4) u (6)3 A a + Ev (7) u (6) 3 A a 

11 1 46 1 1 f 6 

v ± (4) = -0.0148 - j 0.0039 
v ± (7) = 1 

u ± (6) = 0.014 - j 0.027 
X ± = -159.815 + j 18.85 
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b) Eigenvalue s 

A 

X g = —159.815 - j 18.85 (complex conjugate o-f 

A 

c) Eigenvalue X 3 s 

X 3 = -50.0 

AX, = C v (4) u (6) 1 Aa u ^ + C v (7) u,<6)3 Aa 

3 3 3 H 6 3 3 7 6 

v ^ (4) = 0.0003 
v 3 (7) = 0.0039 
u 3 (6) = -0.001 



X 3 = -49.99998 

A 

d) Eigenvalue X H s 



X 

H 


— 


-8.61 + j 


7. 91 


< 

<1 


= 


Cv 4 (4) u h 


(6)] Aa Hfc + Cv s <7> u h <6>] Aa 7fr 


v H (4) 




-0. 1667 - 


j 0. 1322 


v h C7> 


= 


1 




u 4 (6) 


= 


-0.0334 - 


j 0.132 


A 

N 


= 


-8.95 + j 


7.6355 



A 

e> Eigenvalue X s : 

X s = -8.95 - j 7.6355 

A 

i ) Eigenvalue x : 

6 
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< * 



X 

6 


= -1.82 + j 2.02 


AX fr = 


= C v (4) u (6)3 Aa_ + Cv (7) u (6)3 Aa 

<► 6 4 6 6 6 7 6 


< 

* 

X* 


= 0. 1718 - j 0.0161 


v (7) = 


= -0.5737 + j 0.1689 


u fr (6) = 


= 0.01 + j 0.0363 


A 

X 

6 


= -1.846 + j 2.05 


g> Eigenvalue 


A 

X 7 ! 


A 

X 7 


= -1.846 - j 2.05 


h) Eigenvalue 


A 

X 8 ! 


X 8 


= -2.71 


ax 8 - 


= C v (4) u (6) ] Aa_ + Cv (7) u (6)3 Aa 

8 8 46 8 8 76 


v (4) = 
8 


= 0. 1741 


v (7) = 
8 


= -0.6597 


u (6) = 
8 


= -0.0396 


A 

X 8 


= -2.681 



A 

i) Eigenvalue x y : 



S 


= -0.5 


A S ’ 


= C v y (4) u y (6) 3 Aa H6 + Cv ? <7) u ? (6)3 Aa 7fr 


v (4) « 

9 


= 0. 1718 + j 0.016 
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V (7) = -0.5737 - j 0.1689 

9 

u (6) = 0.01 - j 0.0363 

9 

A = -0.4741 

9 

j) Eigenvalue x : 

So =-0-14 

^xo = Cv xo (4) u n (6)1 + Cv io (7) u »o (6,] ^7. 

v (4) = O. 1741 
1 0 

v (7) = -0.6597 
1 o 

ut „ (6 ) = -0.0396 
i o 

X, „ = -o. 1104 

1 o 

For the analysed parameter , our system is 
unsensitive, taking in consideration the small differences 
between the new and original eigenvalues. 

Repeating this procedure for all different parameter 
values will give us the approximation of the new eigenvalues 
and that can be helpfull in a qualitative analysis of what 
parameter has more effect over our model . 

As we want to verify the effect in the time 
response and due to the Software available the analysis will 
be done completely with the calculated exact values. 
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2 . 



Exact Eigenvalues 



Using the Controls Program, we will verify the variation 
of the eigenvalues of the augmented A matrix ( A - FB) with 
respect to the variation of the parameters and by simulation 
we analyse the effects on the time response. 

At the end of the analysis we compare the results of the 
example used for the first order approximation with the 
exact values that were calculated. 

As reference, the eigenvalues of the original system 

ares 



X ± = -159.746 + j 18.979 

X = -159.746 - j 18.979 
2 



X 3 = -49.9992 



X 



H 



X 

X 



5 



6 



X 



7 



8.60697 + 
8.60697 - 
1.81854 + 
1.81854 - 



j 7.90968 
j 7.90968 
j 2.0238 
j 2.0238 



X = -2.70266 
8 

X = -0.503819 
? 

X = -0. 143073 
1 o 
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Each aerodynamic coefficient, that takes part of the 
considered model was changed up to +25 V. and the eigenvalues 

of the perturbed system matrix ( A-augmented ) were analysed 
as well as the TIME RESPONSE of the system . 
a. C 

me* 

On TABLE II we have the eigenvalues related to the 

considered values of C and as we can see, with exception 

m<* 

of the 7— eigenvalue, the changes are small and the system 
is pratically insensitive to that variation. 

The time responses are shown in Figures III. 1 , 111.2,111.3 

and 1 1 1. 4 . 

- Angl e-of -attack (Figure III. 1 > - changes in the same 
direction of C with the minimum value changing about 

moc 

0.5 degrees in the negative direction. The maximum 
value stays almost constant ; 

- Acceleration (Figure III. 2) - the change was less than 

2 

5 ft/sec either the maximum or the minimum values; 

- Velocity (Figure III. 3) - almost the same ; and 

- Altitude (Figure III. 4) - almost the same. 



The change in the eigenvalues are on TABLE III ; the 

time response is on Figure 111.5,111.6,111.7 and III. 8. 

The variation of the eigenvalues are slightly higher 

with changes in C but the system stays stable and the time 

m&p 

response has a small variation. 
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“ Angl e-of -Attack (Figure III. 5) — The minimum values 

increases negativelly about 0.5 degrees but in the 
contrary direction of the coefficient, i.e. , increasing 
the value of C^^ negativelly, the minimum was smaller; 

“ Acceleration (Figure III. 6) - small changes; 

“ Velocity (Figure III. 7) — small changes; 

- Altitude (Figure III. 8) - pratically the same. 



c. 



C 



N<* 



The change in the eigenvalues are on TABLE IV and the 
time response is on Figures 111.9,111.10,111.11 and III. 12. 

The changes in the eigenvalues are small but the effect 
on the time response is more sensible, mainly on the angle- 
of -attack . 

- Angl e-of -Attack (Figure III. 9) - the maximum and the 

minimum increase 1.0 degree in the positive and 
negative directions, respective! ly; 

- Acceleration (Figure III. 10) - maximum changes of 107. 

in the negative direction (minimum) and 207 in the 

positive direction (max i mum) . These values are smaller 

when we increase the value of C fcl : 

Noc 

- Velocity (Figure III. 11) - small changes; 

- Altitude (Figure III. 12) - does not change. 



The change in the eigenvalues are on Table V and the 
time response are on the Figures 111.13,111.14,111.15 and 
III. 16 

- Angl e-of -Attack (Figure III. 13) - Changes almost 1.0 

degree with the maximum and minimum increasing when the 
value of the parameter increases; 
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Acceleration (Figure III. 14) - the maximum increases 

with the value o-f the parameter; 

Velocity (Figure III. 15) — pratically the same -for all 

the considered values o-f the parameter; 

Altitude (Figure III. 16) - constant. 
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TABLE II 

C - SENSITIVITY 

mo: 


C 




-0. 


04 




“0 • 05 








m 
























-159.741 


+ 


J 


18. 9843 


-159.743 




j 


18. 9816 






-159.741 


- 


j 


18.9843 


-159.743 


- 


J 


18. 9816 






-49.9992 








-49.9992 












-8. 55822 


■f 


J 


6.37130 


—8. 56386 


+ 


j 


7. 18736 






-8. 55822 


- 


J 


6.37130 


-8.56386 


- 


J 


7. 18736 






-2.31740 


+ 


J 


2.95574 


-2. 14117 


+ 


j 


2.43275 






-2.31740 


- 


J 


2.95574 


-2. 14117 


- 


J 


2.43275 






-1.79944 








-2. 14200 












-0.51650 








-0.510282 












-0. 14304 








-0. 14306 






I 


c 




-0. 


07 




00 

o 

a 

o 








mot 






















-159.748 


+ 


J 


18. 9762 


-159.751 


+ 


J 


18.9735 






-159.748 


- 


j 


18. 9762 


-159.751 




j 


18.9735 






-49.9992 








-49.9992 












-8.66197 


+ 


J 


8.55389 


-8.72178 


+ 


J 


9. 14021 






-8.66197 


+ 


j 


8.55389 


-8.72178 


+ 


j 


1 

9. 14021 






-1.47981 


+ 


j 


1.81088 


-1.22077 


+ 


J 


1.70081 




-1.47981 


- 


j 


1.81088 


-1.22077 




J 


1.70081 




-3.27131 








-3.67125 












-0.49779 








-0.49148 








i 




-0. 14308 








-0. 14308 
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TABLE III 

C , - SENSITIVITY 
m&p 


C c 

m&p 


-0. 


06 


-0. 


07 








-162.604 


+ j 


16.7131 


-161.200 


+ 


J 


17.9141 




-162.604 


“ j 


16.7131 


-161.200 


- 


J 


17.9141 




-49.9992 






-49.9992 










-5.74603 


+ J 


8.86120 


-7. 13033 


+ 


J 


8.50992 




-5. 74603 


“ J 


8.86120 


-7. 13033 


- 


J 


8.50992 




-1.24694 


+ j 


1.66681 


-1.54230 


+ 


J 


1 . 80335 




-1.24694 


“ J 


1.66681 


-1.54230 


- 


J 


1 . 80335 




—3. 86370 






-3.30548 










-0.49103 






-0.49828 










-0. 14311 






-0. 14309 








C 




—0. < 


09 


-0. 


10 






m&P 
















-158.243 


+ j 


19.9192 


-156.678 


+ 


J 


20.7492 




-15S. 243 


- J 


19. 9192 


-156.678 


- 


J 


20.7492 




-49.9991 






-49.9990 










-10. 1621 


+ j 


6. 98960 


-11.7920 


+ 


j 


5.53519 




-10. 1621 


- J 


6.98960 


-11.7920 


- 


j 


5.53519 




-1.47981 


+ J 


1.81088 


-2.03752 


+ 


j 


2.55411 




-1.47981 


- J 


1.81088 


-2.03752 


- 


J 


2.55411 




-2.27254 






-2.02276 










-0.50786 






-0.51120 










-0. 14307 






-0. 14306 
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TABLE IV 



C N - SENSITIVITY 



* 

z 

u 


0. 


13 


0. 


14 








-159.750 


+ 


J 


18. 9895 


-159.748 


+ 


j 


18.9842 




-159. 750 


- 


J 


18.9895 


-159.748 


- 


J 


18.9842 




-49.9992 








-49.9992 










-8.76667 


+ 


j 


8.27012 


-8.68784 


+ 


J 


8.09494 




-8.76667 


- 


J 


8.27012 


-8.68784 


- 


j 


8.09494 




-1. 18629 


+ 


J 


1.74321 


-1.47163 


+ 


J 


1 . 84903 




-1. 18629 


- 


J 


1.74321 


-1.47163 


- 


J 


1 . 84903 




-3. 63456 








-3.22775 










-0.48878 








-0.49663 










-0. 14314 








-0. 14310 








* 

z 

“J 


0. 


16 


0. 


17 








1 -159.744 


+ 


j 


18.9736 


-159.741 


+ 


j 


18. 9683 


-159.744 




J 


18. 9736 


-159. 741 


- 


j 


18. 9683 


-49.9992 








-49.9992 








-8.52441 


+ 


J 


7.71274 


-8.44077 


+ 


j 


7.50213 


-8.52441 


- 


J 


7.71274 


-8.44077 


- 


j 


7.50213 


-2. 15317 


+ 


J 


2. 32249 


-2.40105 


+ 


J 


2.67918 


-2. 15317 


- 


J 


2.32249 


-2.40105 


- 


j 


2.67918 


i 


-2.20610 








-1.88565 










-0.51041 








-0.51651 










-0. 14305 








-0. 14303 
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TABLE V 

C kl - SENSITIVITY 
N&p 


C M 

N&p 


0.02 


0. 


03 








-159.004 


+ 


j 12.1138 


-159. 748 


+ 


J 


18. 9842 




-159. 004 


- 


j 12.1138 


-159. 748 


- 


j 


18. 9842 




-49.9996 






-49.9992 










-9.37545 


+ 


j 7.11519 


-8.68784 


+ 


j 


8.09494 




-9.37545 


- 


j 7.11519 


-8.68784 


- 


J 


8.09494 




-2.24432 


+ 


j 2.56109 


-1.47163 


+ 


j 


1 . 84903 




-2.24432 


- 


j 2.56109 


-1.47163 


- 


J 


1.84903 




-1.98723 






-3.22775 










-0. 51356 






-0.49663 










-0. 14304 






-0. 14310 








c. . 




0. 


05 


0. 


06 






N oP 














1 


-160. 106 


+ 


j 21.5632 


-160.459 


+ 


J 


23.8437 


I 


-160. 106 


- 


j 21.5632 


-160.459 


- 


J 


23.8437 




-49.9990 






-49.9988 










-8. 25204 


+ 


j 8.25703 


-7.91449 


+ 


J 


8.57438 




-8.25204 


- 


j 8.25703 


-7.91449 


- 


j 


8.57438 




-1.57370 


+ 


j 1 . 85345 


-1.35009 


+ 


J 


1.75196 




-1.57370 




j 1.85345 


-1.35009 


- 


J 


1. 75196 




-3. 18783 






-3.60970 










-0. 49848 






-0.49275 










-0. 14310 






-0. 14312 
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D. CONCLUSION 



As a conclusion for the sensitivity analysis, we have 
the foil owing: 

With respect to the eigenvalues, the pitch control 

system is more sensitive to changes in C and C than 

* m ip N &p 

changes in C and C K1 : and 
m<* Not 

Consi dering the time response, we can verify that the 
angl e-of -attack is sensitive to all the parameters with 
higher variation with C^; the acceleration is more 



sensitive to C and C kl 

mot Not 

value changed less than 
parameters; and the al 

considered var i ati ons. 

It is not necessary to 
analysis for the design of 
only in the eigenvalues 
first order approx i mat i on 



; the velocity has the maximum 

107. with respect to all the 
titude is insensible for the 



emphasize the importance of this 

the control system, but thinking 

and comparing the results of the 

presented for C , with those of 

m&p 



Table III for C changed to -0.06 (AC „ = 0.02) we can have 

m &p m &p 

an idea of the importance of the parameter, using that 
appr ox i mat i on . 
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SENSITIVITY 
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I 
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o 
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Figure III.l Angle-of Attack vs Time 

-for different values of C 
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Figure III. 2 Acceleration vs Time 

for different values of C 
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Figure III. 3 Velocity vs Time 

for different values of C 
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Figure III. 4 Altitude vs Time 

for different values of C 
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Figure III. 5 Angl e-af-Attack vs Time 
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Figure III. 6 Acceleration vs Time 

for different values of C 
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Figure III. 7 Velocity vs Time 

for different values of C 
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Figure III. 8 Altitude vs Time 

for different values of C 
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Figure III. 9 



Angl e-af— Attack vs Time 
-for different values of C 
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Figure III. 10 Acceleration vs Time 

■for different values of C k , 
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Figure III. 12 Altitude vs Time 

for different values 



of 



C 



Noc 



67 



o 




Figure III. 13 Angl e-of-Attack vs Time 

far different values af C^p 
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Figure III. 14 Acceleration vs Time 

for different values of C kl 
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Figure III. 15 Velocity vs Time 

■for different values of C K1 
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Figure III. 16 Altitude vs Time 

-for different values of C„,, 
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IV. ROBUSTNESS ANALYSIS 



A. INTRODUCTION 

In the previous chapter we take into consi der at 1 on the 
change in the eigenvalues and in the time response of the 
system due to variations on the aerodynamic parameters that 
take part into the model; under the robustness analysis we 
verify what we can expect of the system when it is affected 
by any kind of perturbation. 

The main idea is to use the robustness analysis to 
design a new set of feedback gains in order to have the same 
behavior in different situations and/or environments which 
act as perturbat i ons to our missile. 

In classical frequency domain, for single-input single- 
output systems (SISO) , a robust design can be achieved using 
Bode, Nyquist or Nichols plots. With these techniques it is 
possible to define gains that give us gain and phase margins 
for a "robust 1 ' system. For multi-input multi-output systems 
(MIMO) , the classical techniques are no longer valid. 

Taking this in consi derat i on , we have analysed the 
robustness of the original design by means of the minimum 
singular value of the return difference matrix for different 
frequencies and we will use the same technique to improve 
the control system. 
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In order to show why we are applying this technique, we 
start with the Nyquist Criterion for a SISO system and 
extend the analysis to a MIMO system, using the minimum 
singular value theory, as explained by Lehtomaki , Sandell 
and Athans in CRef.61. 

B. THE NYQUIST CRITERION 

Let's take in consideration the SISO system represented 
in Figure IV. 1 . 




Figure IV. 1 Single-Input Single-Output System 

where G(s) includes the controller and plant dynamics 
with unity feedback . 

The Nyquist criterion states that if the open loop 
transfer function G(s) does not have any pole in the right 
half s plane, then the locus of G(s) will not encircle the 
point (—1,0) in the Nyquist plot, where jw is substituted 
for s and the axes represent the real (Re) and imaginary 
(Im) parts of G(jw) for various frequencies. 
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A Nyquist Plot is given in Figure IV. 2 and it can be 



used to design -feedback gains that ensure a robust system. 



The measure o-f the system stability are Gain and Phase 
margin obtained from the Figure, as indicated. 




Figure IV. 2 Nyquist plot - Stable System 

The Phase Margin (PM) is the additional phase lag at the 
gain crossover frequency required to bring the system to the 
limit of instability and the Gain Margin (GM) is the 
reciprocal of the magnitude |G(jw) | at the frequency where 



74 



the phase angle is —180 . The gain crossover -frequency is 
the -frequency at which |G(jw> | is unity. 

The Phase and Gain margins are a measure of how close a 
polar plot is to the -1 + jO point. 

Considering the system subject to additive perturbat i on , 
as we have in the diagram of Figure IV. 3. 




Figure IV. 3 Additive Perturbation 

In order to have the system stable, despite the 
perturbat i on , we need 

|AG(jw) | < II + G ( jw> | (IV. 1) 

As we can see in Figure IV. 2, this condition ensures a 
stable system. 

This idea can be extended to MIMO problems through the 
use of matrix norms and of applying the Multivariable 
Nyquist Theorem. 



75 



C. MULTI-INPUT MULTI-OUTPUT SYSTEM 

The generalization of the SISO theory discussed in the 
previous section has been made for the MIMO problem. 

Let's consider our system as represented in chapter III 
and apply the transf or mat i on in order to have the open loop 
transfer function between input and output. 

x = A x + B u (IV. 2) 

Applying Laplace Transform to the equation (IV. 2), we 

have 

sX = AX + BU (IV. 3) 



where X and U represent the Laplace 

u(t), respect ivel y. Taking the value of X 

X = Csl - AH~ 1 B U 
As Y = C X 

We have Y = C Csl - A] *B 

U 



Transform of and 

in equation (IV. 1) 

(IV. 4) 

(IV. 5) 

(IV. 4) 



This transfer function corresponds to G(s) in Figurt 



IV. 4. 




Figure IV. 4 Return Difference 
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From Figure IV. 4 it is possible to define the "Return 
Difference Matrix" (Kwakernaak and Sivan, [Ref. 63) that 
corresponds to opening the loop at a given point and 
connecting an external input variable at this point; for 
example, at point 1. Assuming the input u as zero, we have: 

Y 4 (s) = - G (s) H (s) V (s) (IV. 5) 

where V (s) corresponds to the new input variable and 
Y i <s) the "returned variable". 

The difference between Y^s) and V(s) is 

V (s) - Y j (s) = C I + G(s)H(s) 3 V(s> (IV. 6) 

The matrix I + G(s)H(s) is defined as the "Returned 
Difference Matrix", with G(s) equal to 

G (s) = CCsI - A3 _1 B (IV. 7) 

As explained by Lehtomaki ,Sandel 1 and Athans in [Ref. 73, 
the Multivariable Nyquist Theorem is derived from the 
relationship 

det [ I + G(s)H(S) 3 = $cl (s ) (IV. 9) 

Sol (s) 

where lei (s) = det (si - A) corresponds to the 

character i st i c polynomial of the open loop transfer function 
and Jol (s) = det (si - A + BF) corresponds to the 

character i sti c polynomial of the closed loop transfer 
function from the system represented in Figure IV. 4. 



77 



The multivariable Nyquist theorem which requires that a 



closed loop stable system have the same number of 
counterclockwise encirclements of the origin by the locus of 
the determinant of I + G(s)H(s) as the number of open loop 
poles that are unstable. 

If I + G<s)H<s) is quasi-singular, a small change in G 
may make the matrix singular ; this causes the detCI + 
G<s)H(s)!3 to become zero and the Nyquist encirclement count 
to change indicating an unstable system. 

Basically, in order to analyse the robustness of the 
system we have to verify how close the return difference 
matrix is to being singular as a function of frequency <o>. 

The natural measurement of the singularity is the 
minimum “singular value", since this is the tightest norm. 

The singular value of a matrix M is defined as 

a. = CX. ( M H M)> 1/2 (IV. 10) 

1 1 

J_| 

M represents the conjugate transpose matrix of a 

l_j 

generic M matrix and any eigenvalue of the product of M 
times M. 

If the minimum singular value is close to zero the 
matrix is quasi -si ngul ar and the system is not robust. 

If one assumes 

*<I + G(jw)H<jw) ) > <IV.11> 
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then it can be shown CRe-f. 7] that the gain and phase 
margins -for the system may be represented by equations 
(IV. 12) and (IV. 13), i. e. , 

Gain margin: GM = 1 (IV. 12) 

1 + 

Phase margin: PM = + cos *11 - (IV. 13) 

2 



Universal 


gain 


and phase 


margins curves, presented by 


Gordon 


CRe-f. 




and repeated 


in Figure IV. 5 allow the 


desi gner 


to 


pick 


a singular 


value that corresponds to a 



desired stability margin -for the system to be designed. 

From singular value plots the we will be able to 
identity the critical -frequencies or range o-f -frequencies 
where the singular value is below the chosen level -from the 
Figure IV. 5. 
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Figure IV. 5 



Universal Gain and Phase Singular Value Plot 
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GAIN MARGIN k (dB) 



D. SINGULAR VALUE ANALYSIS 

Using the program developed by Gordon CRef.8], with 
small changes to make it compatible with the size of our 
system and to input the matrices more easily, we calculate 
the minimum singular value of the return difference matrix, 
with the frequency varying from 0 to 50 rad/sec. 

The program utilizes the IMSL subroutines to generate 
the return difference matrix as well as to made the singular 
value decomposition, given the minimum singular value for 
different f requenci es. 

The return difference matrix depends on the point where 
we open the loop represented in Figure IV. 5. For the 
explained case, considering point 1, we have the "output 11 
difference matrix, that corresponds to I + G(s)H(s); if the 
loop is open at point 2 we have the "input" case and the 
return difference matrix will be I + H(s)G(s). 

From the Figure IV. 5 we choose as reference, for our 
analysis, the minimum singular value as 0.6 that gives us a 

o 

phase margin of about 35 to characterize a "robust" system. 

The minimum input singular value is plotted in Figure 
IV. 6. 

It shows that the system is robust, except -for 
frequencies between 0.6 and 5.0 rad/sec. 
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The results related to the minimum output singular value 
are in Figure IV- 7 and the robustness of the system is poor 
from 0 to 18 rad/sec. 

Those results indicate that the behavior of the system 
under per turbat i ons should present problem for frequencies 
below 10 rad/sec. 

To confirm this expectation, we have plotted the BODE 
diagram (magnitude) for the closed loop system considering 
the angl e-of -attack and altitude with respect to both inputs 
(acceleration and desired altitude). 

The diagrams are on Figures IV-8 to IV. 11, where can see 
the attenuation at the low frequencies indicating that we 
can not expect good response at those frequencies and one 
could further expect problems with robustness as is evident 
from the singular value plots. 

□n the next chapter we try to improve the system 
designing a new set of feedback gains, in order to have a 
robust system- 
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SINGULAR VALUE ANALYSIS 
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Figure IV. 6 Input Minimum Singular Value 
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MIN. SINGULAR VALUE (OUTPUT) 
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SINGULAR VALUE ANALYSIS 
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Figure IV. 7 Output Minimum Singular Value 

l - - - — — 
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Figure IV. 10 Bode Magnitude - Altitude vs. 
Acceleration 
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Figure IV. 11 Bode Magnitude - Altitude vs. 
Desired Altitude 
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FREQUENCY - RADIANS PER SECOND 



V. IMPROVING THE DESIGN 



Under this chapter we try to improve the system by 
calculating new -feedback gains that will yield a more robust 

system, taking into consideration the results indicated in 
Figures IV. 6 and IV. 7. 

The POPLAR program, developed by Gordon CRef.8D will be 
used. A numerical optimization technique is applied in order 
to increase the minimum singular value of the considered 
return difference matrix, therefore, a resultant robustness 
of the design. 

The development in this chapter will be preceded by 
comments on numerical optimization and a description of the 
computer program. 

A. OPTIMIZATION 

The optimization was accomplished by means of the 
Automated Design Synthesis Program (ADS) developed by 
Vanderpl aats CRef . 9D . 

The purpose of ADS as of others numerical optimization 
routines is to find the "best*' possible solution for the 
problem, starting from an initial set of variables and 
updating the design i teracti vel y . The problems can develop 
in convergence of the method and in the computer time 
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needed . 



the 



problem has multiple solutions, 



the 



It 



optimization does not always lead to the absolute optimum. 

The ADS program is designed as a black box optimizer 
which allows the user to choose combinations of one 
dimensional search , optimization algorithm and strategy. 

ADS is used as a subroutine and the parameters that 
correspond to the different applications are chosen by the 
user, as explained in CRef.9J. 

The preferred method used in the ADS applications is 
referred to as Sequential Unconstrained Minimization 
Techniques and can be considered as a method that starts 
with an objective function and the constrains combined into 
an augmented objective function and then minimizing this 
function as if the problem was unconstr ai ned . 

ADS employs penalty function techniques as well as an 
Augmented Lagrange Multiplier. 

As presented in CRef.83, the iteration between the 
user's program and the ADS routine can be represented by the 
block diagram of Figure V.l; where the "info" parameter is 
used as a "flag" for a dialogue between the user and the ADS 
routi nes. 

The details for utilization of ADS are presented by 
Vanderplats in CRef.8D. 
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Figure V. 1 Organization of ADS Program 
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B. POPLAR PROGRAM 



The Gordon's program has the capability of calculating 
the -feedback gains that improve the response o-f the system 
either by pole placement and/or by increasing the minimum 
singular value o-f the return di-f-ference matrix. 

The main objective is to improve the robustness by 
increasing the minimum singular value above a desired level 
chosen -from the universal curve reproduced in Figure IV. 5. 

"Optimum” values of the feedback are calculated by 
minimization of an Objective function whose pole placement 
part is 



OBJ 



V 1 

=Z ( ^ - 



i = 1 



di 






2 

) 



where 

A - real part of the desired eigenvalue i; 

R , . 
di 

X - real part of the computed eigenvalue; 

K . 

1 

A - imaginary part of the desired eigenvalue i; 
ai 

Aj - imaginary part of the computed eigenvalue i; 
i 



and, for the minimum singular value optimization, we 



have 



OBJ t max CO, — jf(jw,p) 

j 



where indicates the desired minimum singular value 

and jj is the minimum singular value at a certain frequency. 
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The optimization procedure changes the feedback gains 
until the minimum singular value is raised to the desired 
1 evel . 

The pole placement and robustness program calculates the 
return difference matrix to be considered, working in the 
complex space (as the analysis is conducted in the frequency 
domain) and the objective functions, calling the ADS routine 
to make the optimization. 

Input's for the program are the matrices of the state 
representation of our missile ( A , B , C and F), where the 
observation matrix (C) was specified as Identity because we 
are assuming all the states are observable as well as 
control 1 abl e. 

Other inputs are the desired minimum singular value, 
desired pole locations and frequency interval to be 
consi dered . 

In our application of POPLAR the initial values of the 
feedback gains are those of the original system and the 
parameters for application of the ADS program are: 

Strategy - Augmented Lagrange multiplier; 

Optimizer - BFGS variable metric method for unconstrai ned 
minimi zation; 

One-Dimensional Search — Polynomial interpolation. 

Two situations will be considered; the INPUT - minimum 
singular value, where the return difference matrix is I + HG 
and OUTPUT-mi ni mum singular value with I + GH as the return 
di f f erence matr i x . 
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C. INPUT - MINIMUM SINGULAR VALUE 

The program was first run to improve the minimum 
singular value for the input situation within the same 
interval of frequencies used in Chapter IV. 

Due to the large CPU time involved, the program was run 
for the critical band of frequencies and the result 
evaluated for the complete interval. 

The minimum singular value taken as reference from the 

o 

Universal curve was 0.6 that implies a phase margin of 35 . 

We have an objective function taking into consideration 
the improvement of the minimum singular value as well as a 
part correspondi ng to the pole placement. The main goal is 
to improve the singular value but with a solution such that 
the poles are in positions where the time response satisfies 
the requi rements. 

The weight of the pole placement part was considered as 
107. of the singular value part. 

The program was started considering only part of the 
feedback matrix as free parameters, increasing this number 
until a reasonable solution was found. 

The best solution was found keeping the original 

feedback gains and calculating gains to feedback all the 
other states to generate the commanded acceleration. 

The computer output of the best solution with the 
corresponding inputs is given in Appendix B. 
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The -feedback gains that improve the robustness of the 



system are the following 



F = 



0.04 -1.15 0.03 -1.7 -9.07 0.32 7.07 0.049 0.986 
000 0000 0 0 



1.4 " 

0 



The plot of the minimum singular value for different 
frequencies is given in the Figure V.2, along with the 
values for the original system. 

As we can see the minimum singular value was increased 
in the critical values of frequency with the minimum 
changing from 0.40 to 0.63 that indicates an improvement in 

o o 

the phase margin from 20 to 35 . Also the bandwidth where 
the system was less robust was changed from 5.0 rad/sec to 
approximately 1.0 rad/sec. 

The robustness of the control system for all the 
considered frequencies is better, arriving close to 0.9 with 
a considerable improvement at low frequencies. 

Using the software Controls, the time response was 
plotted in Figure V.3. 

The maximum and minimum values of the angl e-of -attack 
are smaller than those of the original system; the maximum 
acceleration and velocity are slightly lower and an 
undershoot of about 107. appears on the altitude, that could 
be a problem in some flight conditions. 
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We now 



turn to a consideration of the output return 
difference matrix although this is less important from a 
practical viewpoint. 
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INPUT-MIN. SING. VALUE 
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Figure V.2 - Improvement in the INPUT-Minimum Singular value 
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Improved System 



(INPUT) 



98 




D. OUTPUT 



MINIMUM SINGULAR VALUE 



In the Figure V.4 we have plotted the output minimum 
singular values that correspond to the run where the input 
singular values were improved, as well as the original 
val ues. 

As a considerable degradation o-f the output singular 
values was evident, we try, in the last part o-f this work, 
to arrive to a situation where the system could be 
reasonable robust in both cases, increasing the output 
singular value even with some reduction at the input 
si tuati on . 

A-fter several runs o-f the program, using di-f-ferent 

combinations o-f optimization techniques with di-f-ferent 
starting points, it was veri-fied that due to the 

characteri st i cs o-f the plant, that increases in the input 

singular values resulted in decreases in the output 

situation. 

In order to obtain some improvement, an e-f-fort was made 
to modi-fy the system. 

The -first step was to look at the controllability matrix 
o-f our model , represented by 



1 . 855632D-01 
2.725664D+13 
2.725664D+13 
2. 861 121 D+00 
2. 20181 1D+01 



4.948168D+02 
-2.690449D+13 
-2. 690449D+13 
-2. 822009D+00 
2. 171710D+01 
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-6. 781791D+00 
-2.286975D+00 
9. 855433D+00 
-1 .508382D+01 
4 . 905487D— 0 1 



6. 6891 1 1D+00 
2.255563D+00 
-9. 7 19614D+00 
1 . 487685D+0 1 
-4. 264398D-01 



As we can verify, all the elements are non ero, this 
means that the poles can placed at any desired position; but 
a close look at the 2nd and 3th rows indicates very high 
order numbers compared to the others. This situation 
generates some numerical difficulties in pole placement. 

Using engineering judgement, we have to define what kind 
of changes have to be made for achieving a more robust 
system. 

First we try to get a balanced A-matrix by changing 
units, i. e. , the angles will be in radians instead of 
degrees and the pitch angular rate in radians per second 
instead of degrees per second. Despite some reduction on the 
numbers was not possible to arrive to a satisfactory 
solution . 

More positive effect was obtained by adding the effect 
of the angl e-of -attack in the controller (see Figure II. 3). 

As the original system has the output singular values 
higher than 0.5 for frequencies above 10 rad/sec, we have 
kept the original feedback gains and calculated gains to 
feedback only the angl e-of -attac k , pitch rate and commanded 
actuator . 
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The gains are the -following: 



F = 



0 

0 



0 

0 



0 

0 



0.037 0.027 0 -0.02 0.049 0.986 1.4041 

0 0 0 0 0 0 0 



The result is presented in Figure V.5, compared with the 
original singular values and those -from the improvement in 
the input case. 

The output singular values, compared with the original 
was slightly improved -for -frequencies below 5.0 rad/sec. 

Resultant input singular values are plotted in Figure 
V.6, and they are lower than those -from the improvement in 
the input situation but a little higher than the original. 

The time response were practically unchanged with 
respect to the original model and it is presented in Figure 
V.7. 

Further analysis is needed to obtain a solution with 
high output singular values. 
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OUTPUT - MINIMUM SING. VALUE 




Figure V.4 Output Singular Value Resulted -from the 
Improvement in the Input Singular Value 
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OUTPUT - MINIMUM SING. VALUE 

o 




Figure V.5 Output Singular Value 
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INPUT - MINIMUM SING. VALUE 

o 




Figure V.6 Final Result in the Input Singular Value 
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Final Time Response 
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VI. CONCLUSIONS 



The sensitivity analysis based in the change in 
eigenvalues represents an important role in the design. 

As was demonstrated for the studied missile its behavior 
corresponds to a system insensitive to the variation of the 
considered aerodynamic parameters. 

The analysis gives a good indication of what parameters 
have to be precisely determined. 

The robustness analysis also demonstrated how useful the 
singular value analysis is as an auxiliary tool for the 
desi gner . 

Singular value analysis indicates that the system is 
robust to input perturbat i ons but is deficient in robustness 
to output perturbat i ons. 

The physical nature of the problem indicates that the 
prime concern should be with respect to the input 
per turbat i ons; the output per turbati ons are of minor concern 
but should be kept in mind if unusual conditions should be 
encountered by the missile. 

□n the MIMO design a commomly used method is the Linear 
Quadratic analysis where the performance levels are reached 
by adjusting weighting terms in the "cost 11 function, but the 
results of that method, for non-diagonal R matrices, do not 
necessarily imply a robust system. Using the return 



106 



di f ference matrix we can improve the robustness of the 
design. 

□ne advantage of this procedure is that it permits a 
high level of interaction between the system and the 
desi gner . 

Further problems of a computer nature (i.e., large CPU 
time) are encountered with a high number of states. Analysis 
must be confined to a small frequency range of low singular 
values and the complete system return matrix singular value 
are calculated with the determined feedback gains. 

□ne aspect that should be considered in future 
development is to improve the cpu time used when handle a 
system with a high number of states. 



107 



APPENDIX A 



MISSILE DATA 

Under this Appendix the missile sizing, mass properties 
and aerodynamic parameters are presented as given by Arrow 
[Ref. 11, including the -figures. 

The missile is 1/6 scale o-f the actual circular missile 
configuration and is reproduced in Figure A.l. It is tail 
controlled using four identical control surfaces located 

with +30° dihedral . 

A. GEOMETRY AND MASS PROPERTIES 

In table VI we have the size and mass properties, with 
the respactive values, used for development of the state 
equations. 

Only the uncoupled pitch channel was considered, 
assuming no roll movement. 





Table VI 










MISSILE GEOMETRY 


AND MASS 


PROPERTIES 






Weight 




W 






2525 


lbs 


Length 




1 






168 


in 


Di ameter 




d 






24 


i n 


Re-f erence 


Area -for Coe-f-f. 




S 




IT 


ft 2 


Moment o-f 


Inertia about y 




I 


804 


slug- 


-ft 2 








yy 







108 



\ typical 




Figure A. 1 



Circular Configuration of 



the Missile 
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B. AUTOPILOT /AIRFRAME 

Figure A- 2 shows the block diagram of the autopilot used 
to develop the model of the pitch control channel of the 
mi ssi 1 e. 

The state equations that correspond to the autopilot are 
the following: 

x -150.0 x^ + 150.0 i) 

x = 1.353 x + x - 1.353 x) 

2 13 z c 

x = -6.572 x - 5.0 x + 6.572 f) 

3 1 3 z C 

X = -44.3316 x - 59.11 x 

H 5 6 

x = x - 0.14«2 x, + 0.0395 x, 

5 i 5 p 

x = -188.4 x ^ •+• 188.4 x_ 

x = -0.4608 x - 2.231 x - 0.3406 x + 2.231 x 

7 12 3 H 

- 15.095 x — 20.13 x - 0.1430 x, + 0.4608 d 

S 6 7 ZC 

The input for the autopilot is the commanded 

acceleration i) 

zc 

The values related to the aerodynamic are the following: 

C = - 0.06 
moc 



C 



m&p 



0.08 



C 



No; 



0. 15 



C, 



N&p 



0.04 



1 10 



Dynamic Pressure - q 1650 

Velocity - V 3825.46 



lb/ft 2 
f t/sec 
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tn 
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q(rad/sec) 




APPENDIX B 



POPLAR PROGRAM 



The 

singular 



computer 
value was 



output where the improvement in the input 
found is listed below. 



PILCi LAT 3 PR 0 B OUT A I 



PAoe <oi 



I 



<«« RUN a 5501 >>>>> 



wiax. a. 70000 iNiT preq* a. 70000 step* o.ioooo 



X 10 HT !• 0 . 10000 XIOHT 2 * 0 . 00000 X 1 OMT 3 * 0.00000 



SVMINI* 0.00000 SVM 1 NO* 0.50000 AJ» 0.10000 IDO* 0 
• •• TX A PLANT MATRIX ••• 



- 150.00000 

0.00000 


0.00000 

0.00000 


0.00000 


0.000 00 


■ 705.70003 


0 


ooooo 


0.00000 


1.35310 

0.00000 


o. ooooo 
0.00000 


1.00000 


0.00000 


0 


ooooo 


0.00000 


0 


ooooo 


0.00000 


- 1.57200 

0.00000 


0.00000 

0.00000 


- 5.00000 


0.00000 


0 


ooooo 


0.00000 


0 


ooooo 


0.00000 


0.00000 

0.00000 


0.00000 

0.00000 


0.00000 


0.00000 


-44 


3 3142 


- 57.11000 


0 


ooooo 


0.00000 


0.00000 

0.00000 


0.00000 

0.00000 


0.00000 


1.00000 




14020 


- 0.03750 


0 


ooooo 


0.00000 


0.05000 

0.00000 


o.ooooo 

0.00000 


0.00000 


0.00000 


0 


00000 - 100.37777 


100 


37777 


0.00000 


- 0.44000 

0.00000 


- 2.23000 

0.00000 


- 0.34050 


2.23000 


-15 


07050 


- 20.13000 


“• 


14300 


0.00000 


0.00000 

- 475.00000 


0.00000 

0.00000 


0.00000 


0.00000 


0 


ooooo 


0.00000 


0 


ooooo 


- 50.00000 


0.00000 

0.00000 


0.00000 

1.00000 


0.00000 


0.00000 


0 


ooooo 


0.00000 


0 


ooooo 


0.00000 


0.00000 

0.00000 


0.00000 

0.00000 


0.00000 


0.00000 


-17 


04400 


- 4.70500 


0 


ooooo 


0.00000 



■■■ TX B CONTROL 1XUT MATRIX ••• 



0.00000 0.00000 

- 1.33310 1.53440 

4.57200 - 4.40100 

0.00000 0.00000 
0.00000 0.00000 
0.00000 0.00000 
0.44000 - 0.45450 

0.00000 475.00000 
0.00000 0.00000 
0.00000 0.00000 



■■■ TX C OBSERVATION MATRIX ••• 



1.00000 0.00000 
0.00000 0.00000 
0.00000 1.00000 
0.00000 0.00000 
0.00000 0.00000 
0.00000 0.00000 
0.00000 0.00000 
0.00004 4.00000 

4.44444 4.40044 



0.00000 0.00000 
0.00000 0.00000 
1.00000 0.00000 
0.00000 1.00000 

0.44440 0.40040 



0.00000 0.00000 
0.00000 0.00000 
0.00000 0.00000 
0.90000 0.00000 

1.00444 0.04400 



0.00000 0.00000 
0.00000 0.00000 
0.00000 0.00000 
0.00000 0.00000 

R. 44444 0.44444 
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FJUCi LAT3TR09 OUT 



RAOC »0t 



Al 



0.00000 
0.00000 
0.00000 
a.ooooe 
0.00000 
a. ooooo 
o.ooooe 
a. ooooo 

1.00009 

a.ooooa 

a.ooooe 



0.80000 

9.00000 
0.00000 
0.00000 
a.ooooe 

0.00009 

a. ooooo 
o. ooooo 
a. ooooo 
0.00000 

1.00000 



9.00000 

0.00000 

0.00009 

o.ooeoo 

o.eoooe 



0.00000 

0.00000 

0.00000 

o.eoooo 

0.00009 



0.00000 

0.90000 

0.90009 

0.00000 

o.oeoeo 



1.00000 

9.00000 

o.eeeeo 

0.00000 

0.00000 



o.ooooe 
1 .9000# 
a.aoaoa 

8.00009 

9.90009 



9.00000 

0. 00000 

1. aoooo 
0.00000 
o.eooeo 



'■ twe r reeoeAC* matrix 



0.00000 
o. Tatis 
0.00000 
0.00000 



0.00000 
1 .40440 
0.00000 
0.00000 



0.00000 

o.aiooo 



0.00000 

0.00000 



0.00000 

0.00000 



0.00000 

o.ooooo 



0.00000 

0.00000 



0.04100 

0.00000 



• THE OROEREO COMfLCX EIGENVALUES (INTUT) 



-isi.eoooa 

-151.80000 

-SO. ooooo 
-0.00000 
-0.00000 
-7.00000 

-s. ooooo 
-1.00000 
-1 .00000 
-1 .00000 



-11.09999 

11.01111 

o.ooooo 

-7.10008 

7.90008 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 



•»' THE INITIAL OATA 15 



• TMC 0* OCRED COMMUTED EIGENVALUES 



- 1S9 • 74232 
-151.74232 
-50.02814 
-e.titte 
-B .41448 
-2.44024 
-1.84204 
-1.84204 
-0.50488 
-0.14313 



-10.17508 

18.17508 

0.00000 

-7.81021 

7.81021 

0.00000 

-2.1S070 

2.15078 

0.00000 

0.00000 



THE ORDERED EIGENVECTORS 



-29.44394 


-14.45180 


-29.44394 


14.45188 


1.58392 


8.00008 


-4.43358 


11.39455 


-4.4S3S0 


-11.39455 


•.13420 


0.00000 


4.47925 


-1.32431 


4.47125 


1.32431 


3.40234 


0.00000 


3.79549 


0.00000 










8.24904 


0.09535 


0.24904 


-0.01535 


0.32147 


0.00000 


-1.04880 


-1.34394 


-1.04880 


1.34394 


2.19104 


0.00000 


0.04739 


0.88130 


0.04739 


-0.08130 


0.13223 


0.00000 


0.00472 


0.00000 










-1.30715 


-0.44341 


-1.30715 


0.44341 


-1.54493 


j. ooooo 


5.21377 


4.44100 


5.21377 


-4.44100 


-11.21112 


0.00000 


-4.42382 


-4.25205 


-4.42302 


4.25205 


-8.10271 


0.00000 


-01.51242 


0. 00080 










-0.04045 


-0.34439 


-0.04045 


0.34439 


-0.28442 


0.00000 


-7.73727 


-1.42414 


-7.73727 


1.42414 


1 .37471 


0.00000 


8.72477 


0.53041 


0.72477 


-0.93041 


0.09812 


0.00000 


-8.00145 


o.oooee 










8.00054 


a. 00197 


9.00054 


-0.88197 


8.00551 


0.80008 


0 .54700 


-8.34374 
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FILE i LATSPRD8 OUT 



PACE 003 



*1 



0.54700 


0.34374 


- 0 . 58*03 


0.00000 


- 0.31322 


8.08 299 


- 0.31322 


- 0.08299 


- 0.25444 


0.00000 


- 0 . 2*048 


0.00000 












- 0.01491 


- 1.00104 


- 0.01491 


1.00104 


- 0.24503 


0.00 


0 


- 1.34544 


- 0.94010 


- 1.34544 


0.94010 


0.49479 


0.00000 


0.23844 


- 0.01 


7 


0.23044 


0.01917 


0.19194 


r . ooooo 


0.20211 


0.00000 












- 0.10317 


- 0.15008 


- 0.10317 


0.15008 


- 0.17994 


0.00 


0 


- 1.34340 


- 0.84453 


- 1.34340 


0.84453 


0.49030 


0.00000 


0.23584 




1 


8.23504 


8.02171 


0 . 1914 J 


0.00000 


0.20199 


8.00000 












- 0.04044 


- 0.00944 


- 0.04044 


0.00944 


- 251.15575 


0.00 


0 




• 1.43904 


- 1.28494 


1.43904 


- 12.59149 


0.00000 


2.72214 


0.84 


4 


2.72214 


- 5.84834 


140.74147 


0.00000 


- 1828.24782 


0.00000 












0.00004 


0.00018 


0.00004 


- 0.00018 


0.00042 


0.00 


0 


8.07750 


0.03327 


0.07750 


- 0.03327 


1.32212 


0.00000 


- 0.24593 




1 


- 0.24593 


0.51431 


14.08348 


0.00000 


184.95924 


0.00000 












- 0.00381 


- 0.02884 


- 0.00301 


0.02004 


- 0.02110 


0.00 


0 




- 0.09804 


- 0.40479 


0.89884 


- 3.25274 


0.00000 


- 0.44822 


1.40 


2 




- 1.40442 


- 7.11040 


0.00000 


-24 .47534 


0.00000 













••• THC F FEEDBACK MATRIX 



0.00000 0.00000 

0 . 98 4 1 5 1.40440 

0.00000 0.00000 
0.00000 0.00000 



0.00000 0.00000 

0.01000 0.00000 



0.00000 0.00000 

0.00000 0.00000 



0.00000 0.04900 

0.00000 0.00000 



FREDUENCV 

0.70000 



MIN ADD IN SV 
0.55:84 



2ND ADO IN SV MAX ADD IN SV MIN PROD ADD IN SV 
3.25329 3.25329 1.34110 



FRED SVADMO SVADXD SVH1H SVM1X SVMMO SVMXD 



0.70000 0.05114 35.14094 



0.522202221.19727 2.00000 2.00000 



AAAAA 
A A 



DO DODD SSSSSS 

D D S 

D D S 

D D SSSSS 

D D S 

DO S 

DDDDDD SSSSSS 



FDRTRAN prddram 
FDR 

AUTDHATEO DESIDN SVMTHCSIS 
VERSION 1.00 



CONTROL PARAMETERS 
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met iatsfaub our 



fAOC BOA 



A 1 



1STRAT • J IOFT 

10RAD • 0 NDV 



3 lONEO • I 1FR1NT • loot 

7 NCON • 0 



0FT1H12AT10N RESULT* 



OBJECTIVE FACTION VALIC O.lOt?OC«0* 



design variable* 

LONER 

VARIABLE BOUND 

1 -o.iooooe»o5 

2 -0. 10000E*05 

I -• . 1 0000E *05 

A -0. 10000E*05 

* -0.10000E*05 

* -0.I0000E*05 

7 -0 . 1 0000E*Q5 

FUNCTION EVALUATIONS • 235 

THE OPTIMIZER CONFUTED OUTFUT 15 



THE ORDERED CONFUTED EIGENVALUES 



- 1ST. 7171 1 
-1ST. 71711 

-so.oirss 

-0.04480 
-B.0A4B6 
-7. 0i4*5 
-5.00011 
-0.17A1A 
-0 .17414 
-0 .14105 



-IT. 0*007 
IT. 0*007 
0.00000 
-7.T7487 
7.T7A07 
0.00000 
8.00000 
-0.1T758 
0 . IT 758 
0.00000 



VALUE 

0.1T117E-01 
-0. 114T4E01 
0.30785E-01 
-0.14T51E*01 
-0 . 70 4 T8E*01 
0 . J 1 * 1 4 E • 00 
8. 707ASE*0l 



UFFER 

BOUND 

0. 1 OOOOE*05 
8 . 1 000OE*03 
#. 1 0000E*05 
0.10000E*05 
0. 10000E*05 
i. 10000E *05 
0. 10000E*05 



•••THE ORDERED EIGENVECTORS 



•27.77853 


-12.15010 


-27.77853 


12.15030 


1.71545 


0.00000 


-0.41352 


-0.7*8 42 


-0.41152 


0.T4842 


-2.73255 


0.00000 


-2.22240 


0.00000 


0.44230 


-0.5 7774 


0.442*0 


0.57774 


-3.47572 


0.00000 










0.27374 


0 . 08 28 3 


0.27374 


-0.08283 


0.34835 


0.00000 


0.18800 




0.18800 


0.01511 


14.57874 


0.00000 


-2.87041 


0.00000 


0.03274 


0.00427 


0.012T4 


-0.00427 


-0.00717 


0.00000 










-1.33070 


-0.40257 


-1.33070 


0.40257 


-1.47471 


0.00000 


-0.71773 


0.10034 


-O.T1T73 


-0.18034 


-72.37257 


0.00008 


14.44447 


0.00000 


•0.10777 


0.01774 


-0.I8T7T 


-0.01774 


01.48440 


0.00000 










-0.04024 


-8.S4800 


-0.04024 


0.34800 


-0.30821 


0.00000 


0.50372 


-0.47214 


0 . 503 72 


0.47214 


-1.71182 


0.00000 


-0.80340 


0.00000 


0.02457 


0.00252 


0.0245T 


-0.00252 


•0.00444 


0.00000 










0.00040 


0.00187 


0.00040 


•0.00187 


0.00577 


0.00000 


-0.00022 


0.04230 


-0.00022 


-0.0*230 


0.27422 


0.00000 


0.17743 


0.00000 


•0.032*1 


0.04007 


-0.03241 


-0.04087 


0.24157 


0.00000 










-0.04431 


-0.7*230 


-0.04431 


0.74230 


-0.24534 


0.00000 


0.13517 




0.13517 


0.04574 


-0.43442 


0.00000 


-0.20700 


0.00000 


0.02440 




0.02440 


0.03047 


•0.I7S80 


0.00000 










-0.1 0853 


-0.14151 


-0.10053 


0.14151 


-0.1T40T 


0.00000 


0.12744 




0.12744 


0.04750 


-0.41811 


0.00000 


-0.20337 


0.00000 


0.02440 




0.02440 


0.01044 


-0.17400 


0.00000 










-0.040T5 


-0.00437 


-0.0407* 


0. 00437 


-243.04340 


0.00000 


0.17754 


-0.04027 


0.1T754 


0.04027 


0.77124 


0.00000 


-0.47280 


0.00000 


20.11471 


14.34287 


20.114*1 


•14.34287 


1782.47554 


0.00000 
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P ICE I LATSFAO# OUT 



A 1 



PAOt 001 



0.00005 


0.00017 


0.0000S 


-0.00017 


0.0004* 


0. 00000 


-0.00*13 


0.0040* 


>0.00*13 


-0.0040* 


-0.055 10 


0.00000 


•0.00320 


0.00000 


-2.02*47 


-1.42114 


-2.02**; 


1.42114 


-100.3717# 


0.00000 










-0.00472 


-0.02755 


-0.00472 


0.02755 


-0.02205 


0.00000 


0.01455 


0.01*15 


0.01455 


-0.01(15 


0.31554 


0.00000 


0.422(7 


0.00000 


0.1132* 


1.33711 


0.1132* 


-1.33711 


25.00340 


0.00000 










'■ TX 


F FEEOOACX 


MATRIX ••• 













0.03133 

0.10(15 


-1.14144 0.03070 

1.40440 


- 1 . *15 32 -1.0*174 


9.31531 


7.07442 


0.04109 


0.00000 

0.00000 


0.00000 0. 01000 
0.00000 


0.00000 0.00000 


0.00000 


0.00000 


0.00000 


FREQUENCY 


MIN ADO IN SV 


2 MO ADO IN 3V MAX 


ADO IN SV 


MIN PROD 


ADD IN SV 


0.70000 


0.(7511 


2.9*213 2.0*213 


1.31025 




FREQ 


5VADNO 5VAOXO 


5VM1N 5VM1X 


svmo 


SVMXO 




0.70000 


9.00710 242.01114 


0.710411751.957*7 


2 .00000 


2.00000 
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